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The Baxter limit of adhesive spheres^ is often deemed 
an unphysical model for liquids, due to the thermody- 
namic metastability of its disordered phases with respect 
to crystal formation^. Here we suggest conditions under 
which the disordered phases are actually stable with re- 
spect to the solid for narrow attractive interactions down 
to the Baxter limit. Furthermore, using basic free en- 
ergy considerations we estimate the necessary condition 
for the gas-liquid critical point to be stable with respect 
to crystallization 4 . Possible experimental and simulation 
realizations are briefly discussed. 

For what follows, it is convenient to consider the origi- 
nal Baxter adhesive sphere model, i.e. a square- well fluid 
with well depth e and well width 5, in hard-core diameter 
a units, which we set equal to unity for convenience. Ex- 
cept where explicitly mentioned, the results are for three- 
dimensional systems. For adhesive spheres, the stickiness 
t behaves as a temperature-like quantity, which is defined 
as 



l25e e/kBT , 



(1) 



where k B T is the product of Boltzmann's constant with 
temperature. The equivalent quantity for square- well po- 
tentials with finite width 5 is 



= 4 



((1 + Sf - l) (e £ / fc « T - l) 



(2) 



The second virial coefficient can then be expressed as 
B2v/Bg s = 1 - -L, where B? v s = 2tt/3 is the hard- 
sphere value. 

The critical point of sticky spheres (and, to a good 
approximation, all fluids with a narrow, isotropic attrac- 
tion 6 ) occurs at r c as O.ff 4 . How about freezing? At 
low temperatures, the density of the gas coexisting with 
the solid is very low, p g <C 1, while at high temperatures, 
the freezing density is that of hard spheres. Typically, the 
freezing curve moves rather abruptly, but continuously, 
from one limit to the other. We estimate the temperature 
and corresponding r where this happens and compare it 
with the known r^. 

To estimate the crossover temperature where the den- 
sity of the vapor coexisting with the solid starts to in- 
crease, we make use of the fact that the chemical poten- 
tial per particle /i is the same in the different coexisting 
phases. For the fluid, we use an ideal gas (we check this 
approximation later), while we treat the crystal as a col- 
lection of uncorrelated cells each with free volume Vf. 
Although this last approximation is rather crude com- 
pared to what is usually used in similar systems^, for 
narrow square-well potentials at low temperatures, it is 




FIG. 1: Gas-solid coexistence line from Eq. [4] for square- well 
interaction with Z = 12 and 5 — 0.25, 0.15, and 0.03, from 
top to bottom. Simulation results for the first two cases are 
included for comparison^. 
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FIG. 2: Sketch of the three-dimensional Baxter phase diagram 
following ref.— for the isotropic gas-liquid binodal (dashed 
line) and this work for the gas-solid coexistence curves (solid 
lines). The gas-solid coexistence is found at high r for Z > 6 
and at very low r for Z < 6, while it follows Eq. [7] for Z — 6. 
Gas-crystal coexistence curves for p > 0.1 and gas-liquid co- 
existence for r < 0.1 are drawn for schematic purposes only. 



not unreasonable. This gives 



Ae 
Ze 



As 

k B 



pAv 



2k B T 



k B T 
\np g v f - 1, 



(3) 



where s is the entropy, v is the volume, e is the energy, 
and Z is the coordination number of the crystal. To ob- 
tain the last line, we also assume v to be much smaller 
in the crystal than in the gas, so Av w p~ . Note that 
when A/i — 0, ksT coex /Ze is the only temperature term, 
thus changing the coordination number of the solid ef- 
fectively rescales the coexistence temperature. In other 
words, to a first approximation, the gas coexistence curve 



2 



for a solid with coordination number Zi can be obtained 
from the known gas coexistence curve for a solid with 
coordination number Z\ by multiplying T coex by Z2/Z1. 

Following Sear—, we take vj « (|) for narrow square 
wells in d dimensions. Eq.[3]can then be used to calculate 
gas-solid phase coexistence directly, 
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This is compared to published numerical results of Z = 
12 square- well phase diagrams in Fig. Q] At low p g , to- 
wards the wider limit of narrow wells, there is a fairly 
good agreement between this crude treatment and the 
MC simulated phase diagrams. Results for Z = 4 also 
qualitatively agree with this picture^. 

As a last consistency check, we examine in what regime 



corrections to Av 



, _1 are small. The next order terms 



are the crystal's finite volume, which is of 0(1), and the 
fluid's second- virial coefficient B 2v - Using Eq. 0] as a 
mean- field-like closure for e/T coex gives 



B 2v 
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which for Z > 6 leaves an 0(pg 2 ^ Z ) correction. At low 
densities the two contributing factors are both <C Pg 1 
and have an opposite sign, so they amount to an error in 
T C oex of less than a few percents for narrow well widths 
when p g < 0.1. This level of precision is sufficient for the 
sake of our argument. 

Going to the Baxter limit, 5 — ► 0, we then use the 
definition of stickiness from Eq. [T] to rewrite Eq. |4] as 



1 = - In 



2- d P g5 d - z ' 2 ( 



1 



12t>, 



Z/2 



(0) 



For Z > 2d, S d ^ z / 2 diverges, while at freezing the RHS 
of Eq. [5] should be 0(1). If we fix the density of the gas 
at a small but finite value, it follows that r coex — > 00, 
i.e. the fluid freezes well above r c « 0.11. Conversely, 
when Z < 2d, r coex = 0, so the triple point also lies 
on the t — axis. The interesting case is the marginal 
situation Z = 2d, where we obtain a non-trivial fixed 
crossover T rnPX . For d = 3, the crossover T rnP . r increases 



by less than 15% as 6 varies from 0.25 down to the Baxter 
limit, which is small relative to the divergence obtained 
for Z > 2d. Moreover we expect r coex to lie in the same 
regime as r c . Indeed, in three dimensions Eq. [5] reduces 
to 



(epg) 



1/3 
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(7) 



Using pg = 0.1 as a reference point, we then obtain a 
crossover r coex w 0.03, which is well below r c for the 
isotropic case, as shown in Fig. [2] The gas-liquid critical 
point is thus stable with respect to the crystal. 

Simple cubic crystals, among others, and isostatic 
solids, fulfill the Z — 2d condition for all d. On the 
other hand, monodisperse, isotropically adhesive spheres, 
have invariably Z = 12 crystal phases. However, if ei- 
ther the spheres are not monodisperse or if the inter- 
action is patchy, a lower Z can be obtained. In the 
adhesive limit, even infinitesimally polydisperse solids 
eventually fractionate^, reaching thermodynamic equi- 
librium, which then results in Z rj 12. However, this 
is kinctically nearly impossible and therefore of little ex- 
perimental relevance. Instead, an isostatic polydisperse 
metastable crystal would leave r c stable with respect to 
the gas-solid coexistence. For patchy adhesive spheres Z 
is precisely controlled by the number of patches. Though 
t c also varies with Z^, this effect is unlikely to be suffi- 
cient to bring it below T coex . 

In conclusion, if one manages to reduce the bonding of 
the solid phase down to Z < 2d, the high density fluid 
remains stable with respect to the solid down to the Bax- 
ter limit. Based on rough, but reasonable estimates, the 
three-dimensional gas-liquid critical point would then be 
also stable with respect to crystallization. Furthermore, 
this scheme might be realizable both experimentally and 
through simulation. 
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